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0 Method and apparatus for the acoustic investigation of a casing cemented In a borehole. 



© Method for determining at least one characteris- 
tic of a casing cemented in a borehole, such as 
cement bond and casing thickness, from a reflected 
acoustic signal S(t) obtained by directing an acoustic 
pulse at a substantially nonnal incidence towards a 
radial sector of the wall of tiie casing, said pulse 
stimulating tiiickness resonance within the walls of 
tiie casing. 

The metiiod according to the invention com- 
prises ttie steps of : defining a first time window 

3conresponding to a first portion of signal S(t) includ- 
ing the initial reflection from the casing and subse- 
^quent acoustic returns due to resonance, defining a 
second time window corresponding to a second per- 
^tion of signal S(t) only including tiie initial reflection 
l^from tiie casing, and detennining said casing char- 
ge acteristic from infomnation related to resonance con- 
COtained in said first time window while normalizing 
Qsaid infonmation by Information contained wittiin said 
second time window. 
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METHOD AND APPARATUS FOR THE ACOUSTIC INVESTIGATION OF A CASING CEMENTED IN A 

BOREHOLE 



The invention relates to the acoustic inspection 
of a casing cemented in a borehole which traverses 
undergrourKl strata. It relates more particularly to a 
method for processing a measured acoustic signal 
for fonming signals representative of characteristics 
relating to the casing and to the surrounding ma- 
terials, while simultaneously compensating for vari- 
ations in tiie acoustic pulse generated and the 
acoustic attenuation by the mud filling the 
borehole. 

In general* once a well has reached the desired 
deptii, the borehole is cased, with cement being 
injected into tiie annular space between the casing 
and the wall of the borehole in order to prevent 
hydraulic communication between the various geo- 
logical strata. In order to detenmlne whether such 
unwanted communication nevertiieless exists, mea- 
surements may be performed downhole by means 
of a logging tool in order to determine the quality 
of the bond between die cement and ttie casing. 
Another information of interest which can be de- 
rived from tfiese measurements is the thickness of 
the casing. 

It has long been the practice to use acoustic 
waves for performing such measurements. Such 
techniques have relied on measurements which are 
averages in the circumferential and/or longitudinal 
direction of the casing, and consequently cannot 
identify localized phenomena such as longitudinal 
hydraulic communication paths. A complete pre- 
sentation of prior art techniques can be found in 
U.S. Patent 4;255,798 to Havira, columns 2-7. 

Of ttie techniques ttiat have sought to improve 
the vertical and radial resolution in such inspection, 
tiie technique described in the afc)ove-mentioned 
Havira patent turns out to have been an extremely 
important breaktfirough. This technique consists in 
emitting an acoustic pulse over a radial sector of. 
tiie casing, with the pulse being constituted by 
acoustic waves at frequencies selected to cause 
resonance to appear across the thickness of the 
casing tsetween its outside wall and its in^de wall; 
in determining tiie energy present in a revert)erae 
tion segment of the reflected signal; and in char- 
acterizing from said energy the quality of the bond 
of the cement behind said radial sector of the 
casing. The reverberation segment under consider- 
ation is selected so as to be substantially repre- 
sentative of acoustic revert)eration tietween the 
walls of the casing. Rapid damping of tiie reso- 
nance, in other words low energy, means that ce- 
ment is present behind the casing, whereas slow 
damping, i.e. high energy, means that there is no 
cement 



A logging tool using such a technique Is de- 
scrit)ed in a commercial brochure entitied "Cement 
Evaluation Tool** published by Schlumberger in 
June 1983. The sonde of this tool is centered 

5 Within the casing. 

The sonde includes eight transducers distrib- 
uted helically at 45* intervals, tiiereby obtaining 
good coverage around the periphery of the casing. 
Acoustic pulses are tired sequentially. They are 

10 likewise received sequentially and analyzed and 
transmitted to tiie surface where they are pro- 
cessed. 

In addition, a ninth transducer referred to as 
the "reference" transducer points along tiie axis of 

rs the casing towards a reflecting wall which is plane 
and disposed at a fixed distance from the reference 
transducer. The reflection signal detected by the 
ninth transducer is used to determine in situ the 
propagation time through the borehole fluid (mud), 

20 i.e. the time interval between emission and recep- 
tion of the acoustic wave, and the wave's propaga- 
tion vetocity through tfie mud is deduced there- 
from. Given the propagation velocity of the wave, it 
Is possible to determine the apparent radius of tiie 

25 casing for each of tiie eight transducers. It Is par- 
ticularly advantageous to obtain this radius since It 
makes it possible, in particular, to detect any possi- 
ble deformation of the casing and to monitor the 
centering of ttie sonde inside tiie casing in order to 

30 obtain an indication of tiie validity and the quality of 
tiie measurements as performed. 

In tills technique, a portion is taken from ttie 
reflected acoustic signal, which portion is repre- 
sentative of acoustic reflections t>etween ttie walls 

35 of the casing. 

A signal QC{ is derived therefrom representar 
tive of the quality of tiie bond of the cement with 
tiie casing on ttie basis of ttie energy W2 mea- 
sured in a revert)eration segment of the reflected 

40 signal S. 

The changes in amplitude of such a signal S 
as a function of time and as picked up by a 
transducer are shown in Rgure 2. The effects due 
to the mud are taken into account by normalizing 

45 the measured energy W2 relative to ttie peak am- 
plitude signal W1 of ttie reflected acoustic signal. 
Nevertheless, it turns out that ttiis nonnalization 
does not give complete satisfaction since quantita- 
tive interpretation of tiie measurements shows up 

so problems of divergence in tiie measurements. 

Tests and experiments perfomned by the Ap- 
plicant have shown that ttiese differences can be 
attributed to a large extent to tfie fact that ttie 
conventional processing is sensitive to the prop- 
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erties of the mud, to the characteristics of the 
transducer and of its drive electronics. 

The Applicant has observed that the normaliza- 
tion of ttie portion of tfie signal at energy W2 was 
being performed using a peak amplitude signal W1 
which did not correspond to tfie same frequencies 
as ttiose present in ttie portion of tiie signal cor- 
responding to ttie energy W2, 

In general, ttie energy W1 is not at ttie same 
frequency as ttie energy W2. even in water. The 
energy W1 is a measure of ttie energy maximum 
conveyed by a spectrum component which de- 
pends on the properties of the pulse i.e. the char- 
acteristics of the transducer and its drive eictronics. 
and on the attenuation by the mud. whereas the 
energy W2 includes energy only around Uie reso- 
nance frequency of the casing. 

The Applicant proposes a different concept for 
F>rocessing the reflected acoustic signal, whereby 
variations in acoustic pulse characteristics and miKi 
attenuation can be compensated for in a satisfac- 
tory fashion and ttie reliability of ttie cement bond 
and ottier information obtained from ttie measure- 
ments can be Improved. 

There is provided according to the invention a 
method for determining at least one characteristic 
of a casing cemented in a borehole, such as ce- 
ment bond and casing ttiickness. from a reflected 
acoustic signal S(t) obtained by directing an acous- 
tic pulse at a substantially normal Incidence to- 
wards a radial sector of the wail of the casing, said 
pulse stimulating ttiickness resonance within ttie 
walls of the casing, comprising the steps of : defin- 
ing a first time window conresponding to a first 
portion of signal S(t) including the initial reflection 
from ttie casing and subsequent acoustic returns 
due to resonance, defining a second time window 
conresponding to a second portion of signal S(t) 
only including the initial reflection from the casing, 
and determining said casing characteristic from 
information related to resonance contained in said 
first time window while normalizing said information 
by infbnmation contained within said second time 
window. 

Other characteristics and advantages of the 
invention appear from reading ttie following de- 
scription made witti reference to ttie accompanying 
drawings, in which: 

Rgure 1A is a diagram of an apparatus for 
cement evaluation, tiie apparatus being shown in 
use in a t)orehole; 

Rgure 16 shows in more detail the sonde of 
the apparatus of Rgure 1A; 

Rgure 2 shows a reflected acoustic signal S 
as detected by ttie transducer T, with signal por- 
tions W1. W2 as defined in tiie prior art method; 

Rgure 3 shows a reflected acoustic signal S 
as detected by ttie transducer T. witti signal por- 



tions CW and PW defined in accordance with the 
present invention; 

Rgures 4 and 5 show, respectively, the nor- 
malization window C and the processing window P 
s as defined in the method of the invention; 

Rgure 6 shows ttie frequency spectrum of 
the signal amplitude, with solid lines for the signal 
windowed with ttie processing window and with 
dashed lines for ttie signal windowed with ttie nor- 
10 malization window; 

Rgure 7 shows a portion of the group delay 
spectrum of the signal after normalization, high- 
lighting a dip indicative of a resonance; 

Rgures 8 and 9 show in solid lines the 
15 amplitude spectrum and in dashed lines the group 
delay spectrum of the reflected signal as obtained 
from synttietic data for. respectively, casings of two 
different ttiicknesses; 

Rgure 10 is a diagram of an electromagnetic 
20 model of ttie casing, ttie cement, and ttie mud. for 
use in a mettiod for computing the cement Imped- 
ance and ttiickness of ttie casing. 

With reference to Rgure 1 , an acoustic logging 
apparatus for studying ttie quality of the cement 
2$ ' bond comprises a downhole device 10 suspended 
in a borehole 11 at ttie end of a multiconductor 
cable 12. 

The borehole, filled witti a fluid, typically drill- 
ing mud M. passes through strata 13 and is lined 

30 with casing 14. Cement 15 occupies the annular 
space between the casing and the geological strata 
traversed by ttie borehole. This logging apparatus 
is based on ttie measurement concept of casing 
resonance as disck)sed in the above-mentioned 

J5 Havira patent 

In general, the downhole device 10 is an elon- 
gate body Including a sonde 16, top and tx>ttom 
centralizers 17 and 18 for holding the sonde 16 
centered on the axis of the casing 14. an acoustic 

40 compensation section 22, a hydraulic compensa- 
tion section 32. and an electronics section 19. 

•Various embodiments are possible for the 
sonde 16. In a first embodiment (not shown) the 
sonde may comprise eight acoustic transducers 

45 angulariy distiibuted around the sonde and dis- 
posed helically so as to enable eight sectors of 
casing to be inspected, as described in the afore- 
mentioned Havira patent or U.S. patent 4.685.092 
to Dumont. 

50 In a second embodiment, the sonde 16 has a 
single acoustic transducer T. The sonde 16 has a 
drive motor 20 for rotating tfie assembly carrying 
the transducer T atsout the axis of the sonde. 

Thus, in normal operation, ttie sonde can in- 
55 vestigate ttie entire periphery of the borehole by 
continuously rotating ttie transducer. 

A more detailed description of the structure of 
the sonde in this embodiment is given in US ap- 
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plication N* 494.721 filed March 16, 1990. By 
having a system with an offset drive shaft and an 
abutment it is possible in a calibration mode of 
operation to decouple the reflector 24 from the 
drive shaft so as to cause only the transducer T to 
be rotated Into a position where it looks at the 
reflector. 

The centralizers 18 and 19 are conventional 
and are provided with means (not shown) for estab- 
lishing electrical and/or hydraulic connections be- 
tween the top portion and the bottom portion of 
each of the centralizers. 

The hydraulic compensation section 32 is of 
conventional type and serves to keep the pressure 
of tiie fluid in tfie hydraulic drcuits in tiie downhole 
device 10 at ttie same pressure as ftie borehole 
fluid M. i.e. the hydrostatic pressure of the well, 
thereby avoiding excessive differential pressures 
on sensitive portions of the sonde. 

The section 22 containing ttte acoustic trans- 
ducer has a hollow portion 23 which is open to the 
mud M. The transducer T emits acoustic waves 
through the hollow portion 23 along an axis which 
is substantially perpendicular to the axis of the 
sonde 16. either towards tiie reflector 24 in the 
calibration mode or towards the casing (normal 
measurement mode). The reflector 24 which is 
substantially perpendicular to the emission axis of 
the transducer T simulates a radial sector of the 
casing 14 under Investigation. 

The transducer T is controlled to emit acoustic 
pulses of short duration and having a spectrum 
covering a frequency range from 200 kHz to 700 
kHz. with tiie mean frequency being about 500 
kHz. 

Since ti)e transducer T emits radially, each 
emitted pulse is reflected by tiie various interfaces 
it encounters, namely the casing 14 or tiie reflector 
24. thereby giving rise to an echo signal which is 
detected by the same transducer. 

The detected signals are transmitted to the 
electronics section 19 by appropriate cabling (not 
shown). In the electronics section 19 the signals 
are sampled at a selected rate, put Into digital 
form, multiplexed, and finally transmitted to the 
surface via a modem (not shown) connected to 
conductors in the cable 12. 

On the surface, the conductor's of the cable 12 
are connected to an acquisition system 25 includ- 
ing a modem for decoding tiie transmitted informa- 
tion, and also including a demultiplexer. The sig- 
nals are then stored in a memory 26. The memory 
is connected to a computer 27 suitable for perform- 
ing a sequence of processing stages at each 
depth. These stages are explained in greater detail 
below. 

Signals representative of the depth of tiie de- 
vice 10 are generated by a displacement detector 



comprising a pulse generator 28 controlled by a 
wheel 29 bearing tangentially against tiie cable 12. 
The pulse generator 28 delivers a pulse each time 
tiie cable 12 has moved tiirough a unit lengtii. 

5 The processing performed by the computer 27 
provides a first set of output signals which are 
delivered to a magnetic recorder 30 which stores 
them, for each deptii. in digital form on a magnetic 
tape. The magnetic recording comprises, In par- 

w ticular. the signals delivered by the transducer T 
and it consequentiy enables the data to be pro- 
cessed subsequentiy. Also for each depth, the 
computer 27 provides a second set of output sig- 
nals to an optical recorder 31 which provides a 

75 graphical recording on a film. e.g. of tfie type 
described in French patent 2 491 123. 

Figure 2 shows the variations as a function of 
time In tfie amplitude of a signal S as sensed by 
the transducer T. and illustrates tiie metiiod of 

20 windowing signal S disclosed in tiie Havira patent 
A signal representative of tiie quality of tiie 
bond between the cement and tiie casing is pro- 
duced from the energy W2 measured in a portion 
of tfie reflected signal S which essentially cor- 

25 respond to revert^erations between tiie walls of the 
casing as excited by the pulses emitted by tiie 
transducer. 

In tiie window ta - ta, ttie energy W2 may also 
include energy due to multiple reflections on tiie 

30 geological strata. These reflections can be ot>- 
served from instant ti . 

The energy W2 is nonmallzed by the maximum 
peak energy W1 in tiie sensed signal S. This 
normalization serves to eliminate effects due to the 

35 mud, but it Is not fully satisfactory as explained 
hereinabove. 

Rgure 3 shows the same variations in am- 
plitude of the signal S as a function of time and as 
sensed by the transducer T, togetiier with two 

40 windows defined as explained below, whtoh both 
include a time reference or time origin to located In 
a portion of tiie signal which corresponds to tiie 
Initial reflection from the casing wall. 

45 

Determining a time origin 

According to tiie invention, a stable time origin 
is determined so as to escape botii from noise and 
50 from variations in transducer response over dif- 
ferent measurements. 

To do this, several methods can be contem- 
plated. 

In a first possibility, an analytic signal Sa(t) is 
55 determined as foltows: 
Sa(t) = S(t)-IH(S(t)) 

where S(t) is the measured signal and H is the 
Hilbert transform of tiie measured signal, with i 



4 



7 



EP 0 395 499 A2 



8 



being the square root of -1. The Hilbert transform is 
obtained by any known algorithm. 

The time origin to is defined by the instant at 
which the analytic signal has a peak. 

Another method of obtaining the analytic signal 
Sa(t) consists in performing a Fourier transform by 
any fast Fourier transfonm algorithm, then in remov- 
ing the negative frequencies by filtering, and sul>- 
sequentiy in performing a further fast Fourier trans- 
form on the signal obtained after filtering. 

The instant at which a peak occurs in the 
analytic signal is taken as being the origin to. 

An even faster method consists in rectifying 
the signal received by the transducer, in extracting 
the amplitude maximum from the rectified signal, 
together with two extreme values on either side of 
said maximum so as to obtain three points defined 
by (time, amplitude) coordinates. 

A second order polynomial passing through 
these three points is then determined in conver)- 
tional manner. 

The maximum of the resulting polynomial is 
determined and the time coordinate to of ttie maxi- 
mum is taken as being the time origin and serves 
as a reference for determining the position and 
duration of the windows. 



Windowing the signal S(t) 

A first portion PW of the signal S(t) is defined 
by a wide time window (-ti, ta) designated as 
processing window hereinafter. The signal portion 
PW obtained in this way includes the initial reflec- 
tion of the acoustic signal from Jtiie inner wall of the 
casing and acoustic returns due to revert>erations 
between the walls of said casing caused by its 
resonance in response to tfie acoustic pulse. 

A second portion CW of signal S(t) used for 
normalization purposes Is defined by a narrow time 
window (-ti. +ti). also called the normalization 
window, which window is centered at time origin to 
and can be defined as the early portion of tiie wide 
window. The signal portion CW thus defined mainly 
comprises ttie initial reflection from tiie wall of the 
casing. 

Rgures 4 and 5 respectively show an example 
of the windowing signal C(t) corresponding to a 
narrow window as used for normalization, and an 
example of the windowing signal P(t) correspond- 
ing to a wide window as used for perfonnning tiie 
processing descrit>ed below. 

Windowed signals corresponding to the por- 
tions CW and PW may be obtained by multiplying 
tiie signal S(t) by tiie signals shown in Figures 4 
and 5, namely tiie narrow window signal C(t) and 
the wide window signal P(t). 

The time intervals spanned by tiie windows are 



defined as multiples of the nominal resonance pe- 
riod of the casing, so as to automatically adapt the 
windows to the type of casing under investigation. 
The time intervals -ti.to and to.ti should lie in 

5 the range t)etween 1.5 and 3 times the nominal 
resonance period of tiie casing, and are preferably 
equal to 2.5 timois this nominal period. The nor- 
malization window -ti .tt tiius lasts for a period lying 
between 3 and 6 times tiie nominal resonance 

70 period of tiie casing, preferably a period equal to 5 
times this nominal period. 

A suitable range for tiie time interval tb.ts of the 
processing window is from 6 to 12 times, prefer- 
ably 6 to 8 times tiie nominal resonance period of 

75 the casing. Thus« the entire time Interval -ti.t^ 
spanned by the processing window is suitably be- 
tween 7.5 times and 15 times, preferably between 
7.5 and 1 1 times the nominal resonance period of 
tiie casing. The nominal resonance period of the 

20 casing is determined from tiie nominal tiiickness of 
tiie casing and tiie sound velocity in tiie casing, 
botii parameters being known. 

For tiie window signals C(t) and P(t), Hamming 
windows can be used. Figure 4 shows such a 

25 Hamming window for C(t). centered on time origin 
to. As regards the processing window P(t), although 
Rgure 5 shows for it a symmetrical curve, it is 
preferred to use a Hamming window centered on 
tiie time origin tb but truncated at time ti i.e. an 

00 asymmetrical window signal. 



Characterizing ttie resonance of ttje casing 

35 The resonance of the casing in response to the 
acoustic pulse depends on the cement-casing 
bond, which will be expressed as an acoustic im- 
pedance of the cement and tiie casing tiiickness, 
which are the primary parameters of interest In the 

40 casing inspection. These parameters can be deter- 
mined from a characterization of a resonance mode 
of the casing e.g. tiie fundamental mode. 

As noted above, the wide processing window 
P(t) is defined so as to contain the direct reflection 

45 from the casing and the early part of tiie reso- 
nance. The narrow window C(t) only contains the 
reflection from the casing and as such, tfie win- 
dowed signal CW(t) provides information on the 
system formed by the pulse generating assembly 

so transducer and drive electronics - and the mud. In 
ottier words, it provides the nssponse of tiie 
"infinite-block" {i.e. a casing of infinite thickness 
and extent). In tiie determination of tiie parameters 
of interest, the information contained in the signal 

55 portion CW will be used to nonmalize the informa- 
tion contained in signal S(t) within tiie processing 
window P(t) so as to automatically compensate for 
variations (with temperature, time etc) in the char- 
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acteristics of the acoustic pulse and tiia attenuation 
of the mud. 

In order to find and characterize the resonance, 
the frequency spectra of the windowed signals PW 
and CW are used. Rgure 6 shows in solid lines the 
amplitude spectrum PV^io) obtained from signal 
PW(t). using a discrete Fourier transform, and in 
dashed lines ttie amplitude spectrum CW(») ob- 
tained likewise from signal GW(t). The resonance of 
the casing appears as a marked dip in signal PW- 
(«). 

Now. the amplitude spectrum being highly de- 
pendent on transducer variations and attenuation 
by tiie tx>rehole fluid, it is preferred to use tiie 
group delay spectrum. The group delay is defined 
as the derivative of radian phase vs. angular fre- 
quency and is computed from ttie phase spectrum. 
The group delay response is advantageous in that 
it tends to t>e smooth and linear except at resonan- 
ces and ttierefore. allows easy Identification and 
characterization of the resonance. 

This is illustrated in Rgure 8. which shows the 
theoretical amplitude (solid line) and group delay 
(dashed line) responses of a casing of 178 millime- 
ters (7 inches) diameter and 4 millimeters thick- 
ness, in the case of a solid cement bond (cement 
impedance of 8 MRayi). Rgure 9 shows likewise 
the amplitude and group delay responses of a 
casing of 178 millimeters diameter and 16 mm 
thickness, witti tiie same cement bond. 

Rgure 7 highlights tiie group delay spectrum 
around a resonance. The resonance can be char- 
acterized by the central frequency Fo of the reso- 
nance, corresponding to a minimum group delay 
rmin. the depth H of ttie dip conresponding to tiie 
resonance, an* the fractional bandwidth B/Fo. The 
t>andwidth B is tiie width of the dip at a level Ar 
above the group delay minimum rmm equal to a 
predetermined percentage. e.g. 40%. of the depth 
H of the dip : B - F2 - Fi. with the corresponding 
group delays t2 and n such that ra - Train = ri - 

rmin = At = 0.4H. 

The central frequency F© Is primarily depen- 
dent on ttie casing thickness while the depth and 
fractional bandwidtti are primarily influenced by the 
acoustic impedances of tiie mud and the cement 



Determining casing thfckness and cement imped- 
ance 

A normalization step of the group delay spec- 
trum is first performed in order to eliminate any 
slope, by linearizing tiie group delay spectrum of 
the normalization signal CW and subtracting it from 
the spectrum of ttie signal PW. 

Then the parameters of the resonance are de- 
termined, by searching a minimum in tiie normal- 



ized group delay spectrum within a pre-determined 
frequency range kx:ated around the nominal reso- 
nance frequency of the casing. This yields values 
for ttie parameters FO. H and B/FO. which will be 

5 used as tiie "measured parameters'* in ttie follow- 
ing description. 

After the measured parameters of the reso- 
nance have been determined from the normalized 
group delay spectrum, tiie casing ttitckness and 

to cement impedance are determined by iteration us- 
ing a simple model, shown in Rgure 10. This 
model is a planar model (i.e. ttie casing wall is 
modelled as a plane) defined by the cement im- 
pedance Zee and ttie casing ttiickness d. with mud 

75 impedance and casing impedance Zp. Mud 
impedance Zm is obtained by separate calibration 
and casing impedance Zp is a known data Resis- 
tances R can be added to tiie model to perform 
corrections taking into account the curvature of ttie 

20 casing, the resistances being set to values defined 
in accordance with the diameter of the casing. 

The response of the model to a Dirac pulse e 
is detenmined in the frequency domain, with the 
initial value for tiie cement impedance Zee obtained 

25 from the measured depth H of resonance by using 
look-up tables, and tiie initial value of casing ttiick- 
ness derived from central frequency Fo. This yields 
an impulse response spectrum of the model. 

This spectrum is ttien multiplied by tiie nor- 

30 malization spectrum CW(«) representative of ttie 
transducer-mud system to obtain what can be 
tenmed the transducer-mud spectrum. 

The processing window Is then applied in the 
frequency domain by convolution with the 

35 transducer-mud spectrum, and the group delay 
spectrum of the convoluted signal is computed. 

This group delay spectrum is again normalized 
as described at)ove by subtracting tiie linearized 
group delay spectrum of the normalization signal 

40 from the group delay spectrum of the convoluted 
signal. 

The normalized group delay spectrum is then 
used to detenmine "model" values for the param- 
eters of the resonance, namely the central fre- 

45 quency Fmod corresponding to the group delay 
minimum t„m identified in tills spectrum within ttie 
specified frequency range, and tiie fractional band- 
width Bmod taken at a level rmod = tmid -f* Ar (Ar 
being as defined at)ove). These model parameters 

50 are compared with the measured parameters. From 
the comparison, new estimates for the cement im- 
pedance and casing tiiickness are defined for a 
second iteration. The iteration is continued until ttie 
model group delay response is characterized by 

55 the same parameters as the measured resonance. 
Iteration may be performed by Newton-Raphson 
fitting technique, which typically requires only three 
iterations. 
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The planar model used in this method is quite 
simple and has the advantage of minimizing pro- 
cessing time. 

A more exact model, based on the cylindrical 
geometry of the casing, could be used instead of 
the planar model. Such a cylindrical model would 
include additional parameters (casing diameter, 
transducer dimensions, transducer offset from the 
borehole axis, sound velodty in mud). 



Alternative computation methods 

Other computation methods can be envisioned 
with the same basic definitions of time windows C- 
(t) and P(t) as specified above. An alternative meth- 
od consists of the following steps. 

The entire signal S(t) is converted into the 
frequency domain before any windowing. The 
spectrum S(«) thus ot)tained is normalized by the 
spectrum CW(«) derived from the calibration signal 
CW(t). the normalized amplitude being the ratio of 
the amplitudes and the normalized phase the dif- 
ference between the phases. 

The processing window is then applied to the 
normalized spectrum In the frequerrcy domain, by 
convolution of the normalized spectrum with the 
spectrum of the processing window. 

The group delay spectrum is then calculated 
fifom the result of the convolution. This allows the 
resonance to be characterized as described above, 
thus yielding central frequency, depth and frac- 
tional bandwidth. 

From this point, it is possible to determine the 
cement impedance and casing thickness by direct 
computation, using pre-established look-up tables 
relating these parameters with the parameters 
which characterize the resonance. Another ap- 
proach Is to use the model as defined above. The 
computation steps will be as explained above, ex- 
cept that the step of normalizing the group delay 
spectrum can be dispensed with. 



Claims 

1. A method for determining at least one char- 
acteristic of a casing cemented in a t)orehole, such 
as cement bond and casing thickness, from a re- 
flected acoustic signal S(t) obtained by directing an 
acoustic pulse at a substantially normal Incidence 
towards a radial sector of the wall of the casing, 
said pulse stimulating thickness resonance within 
the walls of the casing, comprising the steps of : 
defining a first time window corresponding to a first 
portion of signal S(t) including tiie initial reflection 
from the casing and subsequent acoustic returns 
due to resonance, defining a second time window 



corresponding to a second portion of signal S(t) 
only including ttie initial reflection from the casing, 
and detenmining said casing characteristic from 
information related to resonance contained in said 
5 first time window while normalizing said information 
by information contained within said second time 
window. 

2. A method according to claim 1 . wherein said 
second time window is defined so as to include the 

10 amplitude maximum of signal S(t). 

3. A method according to claim 2. wherein the 
definition of tiie first and second time windows 
comprises the steps of determining the instant at 
which the amplitude maximum occurs in the signal. 

/5 said instant being taken as a time origin tb. and 
defining first and second time windows on the 
basis of this time origin. 

4. A metiiod according to claim 3, wherein, for 
determining the time origin, an analytic signal Sa(t) 

20 of the response S(t) is determined: ttie amplitude 
maximum of tfie analytic signal is determined; and 
tiie time origin to is determined as being tiiat in- 
stant at which the amplitude maximum occurs in 
the analytic signal Sa(t) 

25 5. A method according to claim 4. wherein the 
analytic signal Sa(t) is obtained by a Hilbert trans- 
form. 

6. A method according to claim 4. wherein the 
analytic signal Sa(t) is obtained by performing a 

30 Fourier transform, filtering to eliminate negative fre- 
quencies from the spectrum, and performing a 
further Fourier transform on tiie signal as obtained 
after filtering. 

7. A method according to claim 3. wherein the 
35 time origin is defined by performing the following 

steps: extracting the absolute values of the am- 
plitude maximum in the acoustic signal and of two 
extreme values situated on either side of said maxi- 
mum; determining a second degree polynomial fit- 
40 ting these three values; and defining the time origin 
as tiie instant at which the polynomial is at a 
maximum. 

8. A metlKxl according to claim 3. wherein the 
time intervals spanned by the time windows are 

45 defined by multiples of the nominal resonant period 
of the casing. 

9. A metiiod according to claim 8. wherein the 
second time window is centered on tiie time origin 
to, and extends on eitiier side of tiie time origin 

so over a period of time lying in the range from 1.5 
and 3 times tiie nominal resonance period of ttie 
casing. 

10. A method according to claim 8 or 9. 
wherein tfie first time window extends Ijefore the 

55 time origin to over 1.5 to 3 times tiie nominal 
resonance period of the casing and after tiie time 
origin over 6 to 12 times said nominal resonance 
period. 
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11. A method according to claim 1, wherein 
said characteristic of the casing is determined by 
combining infomnation derived from the frequency 
spectrum of the acoustic signal S(t) within the first 
time window and information derived from frequen- 5 
cy spectrum of the acoustic signal S(t) within the 
second time window. 

12. A method according to claim 1 1 comprising 
the step of determining the frequency spectra of 

the first portion of signal S(t) and of the second * lo 
portion of signal S(t). respectively, normalizing the 
spectrum of the first portion by the spectrum of the 
second portion, and characterizing from said nor- 
malized spectrum the resonance of the casing, the 
parameters of the resonance being indicative of the is 
characteristics of the casing. 

13. A method according to claim 12, wherein 
the frequency spectra are the group delay fre- 
quency spectra. 

14. A method for acoustic inspection of a cas- 20 
tng cemented in a borehole, comprising the steps 

of obtaining a reflected acoustic signal S(t) by 
directing an coustic pulse at a substantially normal 
incidence towards a radial sector of the wall of the 
casing, said pulse stimulating thickness resonance * 2s 
within the walls of the casing, defining a first time 
window conresponding to a first portion of signal S- 
(t) including the initial reflection from the casing 
and sut)sequent acoustic retums due to resonance, 
defining a second time window corresponding to a 30 
second portion of signal S(t) only including the 
initial reflection from the casing, and detenmining at 
least one characteristic related to the casing from 
information related to resonance contained in said 
first time window while nonmalizing said information 3s 
by information contained within said second time 
window. 

15. An apparatus for acoustic inspection of a 
casing cemented in a borehole, comprising a 
sonde adapted for displacement in the borehole 40 
which includes transducer means for obtaining a 
reflected acoustic signal S(t) by directing an acous- 
tic pulse at a substantially normal incidence to- 
wards a radial sector of the wail of the casing, said 
pulse stimulating thickness resonance within tiie 45 
walls of the casing; and signal processing means 

for performing the steps of defining a first time 
window corresponding to a first portion of signal S- 
(t) including the initial reflection from tiie casing 
arfd sut>sequent acoustic retums due to resonance, so 
defining a secorKi time window corresponding to a 
second portion of signal S{t) only including the 
initial reflection from the casing, and determining at 
least one characteristic related to the casing from 
information related to resonance contained in said 55 
first time window while normalizing said information 
by Information contained within said second time 
window. 
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